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STUDY OF CHUGGING INSTABILITY WITH LIQUID-OXYGEN 
AND GASEOUS-HYDROGEN COMBUSTORS 
by Marcus F. Heidmann, Daniel E. Sokolowskl, and Lawrence A. Diehl 

Lewis Research Center 

SUMMARY 

Chugging instability boundaries for a small-scale, liquid-oxygen and gaseous- 
hydrogen combustion system were experimentally determined over a broad range of oper- 
ating conditions to establish quantitatively the experimental effect on stability of 
combustion- system variables. Parallel-jet injection of oxygen into diffusely injected hy- 
drogen was used. Experimental variables included the characteristic combustor length 
and, for the oxygen system, the number of jets, injector orifice length to diameter 
ratio, flow- system volume, manifold pressure loss, and atomization mechanism. The 
experimental results were compared by using the oxygen flow- system pressure drop and 
the range of stable operation as criteria. 

The experimental results were evaluated with the aid of an analytical model that 
characterizes the combustion process by a time lag (delay time between injection and 
burning) and an interaction index (sensitivity of combustion processes to pressure oscil- 
lations). The analysis showed the stability behavior of the combustor system to be rel- 
atively insensitive to the response time of the oxygen flow system. Gas entrainment in 
liquid oxygen due to local cavitation and/or heat transfer was identified. The average 
interaction index for all tests was in the range 0 to 0. 2. A qualitative model for the 
breakup time of liquid jets was used to evaluate time-lag variations. Time lag and jet- 
breakup time appeared to be related to (1) the jet-formation process, which varied with 
configurations, and (2) the flow properties of the liquid jet and surrounding gas, which 
varied with operating conditions. An empirical equation based on typical jet- breakup be- 
havior was formulated to give time-lag variations over a range of combustor conditions. 


INTRODUCTION 

The results of an experimental study of low-frequency (chugging) instability in a 
nominal 200-pound-thrust liquid-oxygen and gaseous-hydrogen combustion system are 



presented. The study was initiated to investigate flow- system and combustor variables 
that were known to have or suspected of having an effect on the stability of small-scale 
experimental combustor systems. Variables in addition to propellant mixture ratio and 
total flow rate include characteristic combustor length, injector orifice area, injector 
orifice length to diameter ratio, injector manifold volume, supply line length, and mani- 
fold pressure losses. For most of these variables, an effect on stability had been noted, 
although not fully documented, in previous small-scale experimental programs. In some 
instances, the effect on stability appeared unexplained or contrary to previously estab- 
lished principles for stability. These variables specifically emphasize the liquid-oxygen 
flow system. The cryogenic flow system is susceptible to local cavitation and may cause 
abnormalities in stability characteristics. The cryogenic flow system has been largely 
ignored in previous experimental studies of chugging instability. A secondary objective 
of this study was to determine the probable deviations from pure liquid flow. 

In the present study, experimental data were acquired to establish quantitatively the 
effect of the previously mentioned variables on stability. Changes in stability are com- 
pared with changes in experimental performance parameters. An analytical model rep- 
resentative of the combustor system is also developed, and the characteristic stability 
parameters provided by the analysis are used to evaluate the experimental results. 
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SYMBOLS 

o 

effective area of flow- system supply line, in. 

2 

resistance parameter of flow system, (C^A^/C^A^ 

o 

area of injector orifice, in. 

2 

area of upstream orifice, in. 

flow coefficient of injector orifice, m^/A. (2gp AP.)^^ 
flow coefficient of upstream orifice, m ( - ) /A u (2gp AP.)^ 
characteristic exhaust velocity, ft/sec 
elasticity parameter of flow system, (2 AP.pZAe)/mQ0g 
frequency, cps 

2 

gravitational constant, ft/sec 

inertia parameter of flow system, Zm 0 /(2 AP i gA0 g ) 
flow resistance, AP t /(m 0 ) 2 , sec 2 /(in. 2 )(lb) 
characteristic length of combustor, in. 
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length to diameter ratio of injector orifice 
function notations 

length of flow- system supply line, in. 
mass flow rate, lb/sec 
mass flow rate of hydrogen, lb/sec 
mass flow rate of oxygen, lb/sec 
interaction index 

oxidant to fuel mass flow rate ratio 

combustor pressure, psia 

hydrogen injector orifice pressure drop, psi 

oxygen injector orifice pressure drop, psi 

total oxygen-flow- system pressure drop, psi 

pressure-drop parameter of injector orifice, P c /(2 AP^) 

universal gas constant, 1545 ft-lb/°R(lb mole) 

Laplace transformation parameter 
combustion gas temperature, °R 
time, sec 

velocity of unburnt hydrogen in combustor, ft/sec 
oxygen j et velocity, ft/sec 

position of concentrated flow- system elasticity 
dimensionless time, t/0 

6 

dimensionless time lag, r/0 

_ i 

compressibility of flow- system fluid, psi 
characteristic velocity efficiency 
mean gas residence time in combustor, sec 
real part of Laplacian operator when s = X + io> 
fractional flow perturbation, j(mQ - (z) 

fractional flow perturbation with retarded action, 

O 

density of oxygen, lb/ft 
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o 

p density of combustion gases, lb/ft 

D 

r total time lag from instant of injection to instant of combustion of a propellant ele- 
ment, sec 

Tj liquid-jet-breakup time, sec 

cp fractional pressure perturbation, 

cp' fractional pressure perturbation 

oo dimensionless frequency 

Superscript: 

(~~) mean steady- state value 


[( P c - P c)/ P c>> 

ith retarded action, j(P c - P c )/pJ(z - 6) 


COMBUSTION AND FLOW SYSTEMS 
Test Facility 

Various combustor configurations were evaluated in an experimental rocket test fa- 
cility equipped with gaseous-hydrogen and liquid-oxygen flow systems. A schematic dia- 
gram of these flow systems is presented in figure 1. Both the fuel and oxidizer systems 
utilized a pressure-regulated feed to control flow rates. The mass flow rate of hydrogen 
to the combustor was metered by a critical-flow venturi. Liquid oxygen was supplied 
from a pressurized tank immersed in a liquid-nitrogen bath. A cooling coil, turbine- 
type flowmeter, and propellant valve were also immersed in liquid nitrogen to assure a 
constant low-temperature supply of oxygen. Nominal flow rates ranged as high as 
1. 0 pound per second for oxygen and 0. 25 pound per second for hydrogen. Lengths and 
volumes of feed lines are also shown in figure 1. The length of the liquid-oxygen supply 
line past station 3, hereinafter referred to as the liquid- oxygen manifold, is discussed in 
greater detail in a later section. 


Combustor Configuration 

The experimental rocket combustor, shown in figure 2, was nominally rated at 
200 pounds thrust and a chamber pressure of 400 pounds per square inch absolute. An 
uncooled copper combustion chamber, with an internal diameter of 2 inches, was fabri- 
cated in sections to facilitate interchanging several chamber lengths. The nozzle was 
constructed with a 3-inch convergent section and a water-cooled throat. A mean nozzle- 
throat diameter of 0. 62 inch gave a contraction ratio of 10. 4. As shown in figure 3, 


4 



gaseous hydrogen was injected in a diffuse manner through a multihole plate that spread 
the hydrogen throughout the chamber cross section. This hydrogen system was held con- 
stant, and design parameters for the liquid- oxygen system were changed. The liquid- 
oxygen manifold and injector are shown in figure 4. A variable number of closely spaced 
rows of 0. 040- inch- diameter jets consisting of 20 jets per row produced a showerhead 
stream parallel to the longitudinal axis of the chamber. The resultant pattern was essen- 
tially a sheet of liquid oxygen surrounded by gaseous hydrogen. With this design, the 
manifold and orifice configuration could be changed without substantially affecting the 
basic injection pattern or chamber geometry. 


Reference Configuration 

The configuration changes selected for study were deviations from a basic or refer- 
ence configuration. The chamber length for this reference configuration was chosen to 
given an L* of 156 inches, while the liquid- oxygen injector and manifold are as depicted 
in figures 5 and 6(a). The basic injector consisted of two rows of 20 equally spaced ori- 
fices. All design parameters that were studied were variations of the reference config- 
uration. 


Combustor Configuration Variables 

Characteristic length of combustor . - Characteristic length L* was varied to study 
changes in the mean gas residence time 9 . The residence time is related to L* by the 
relation 9g = L*c*/gRT • This residence time is a characteristic parameter of an un- 
stable combustor (ref. 1). Characteristic combustor lengths of 81 and 224 inches were 
used to study the effect of both increasing and decreasing 6 g. 

Length to diameter rati o of injector orifice . - Jet-length variations (figs. 6(b) and (c)) 
were examined for changes in orifice flow characteristics. Two longer L/D ratios of 
3. 90 and 6.25 were compared with the reference L/D of 1. 56. 

Number of jets . - The number of oxygen orifices (jets) was increased to 60 and de- 
creased to 20 to vary the distribution of pressure drop and jet velocity at a constant total 
flow rate. These changes are depicted in figures 6(d) and (e). 

Flo w- system volume . - The volume of the liquid- oxygen cavity (fig. 6(f)) was in- 
creased from 0.25 to 1. 56 cubic inches by increasing the cavity diameter from 1/4 to 
5/8 inch. The larger cavity increased the flow- system capacitance at the injector and 
also gave lower values of jet- entrance cross velocity; both these could give significant 
effects on chugging stability. 
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The volume of the supply manifold was increased from 2. 85 to 3. 90 cubic inches by 
extending the manifold side branch from 5 to 12 inches. This increased volume gave a 
significant change in overall flow- system capacitance that could affect stability. 

Diameter of upstream orifice. - Upstream orifice diameter was reduced from 0. 070 
to 0. 052 inch to test the effect of a change in pressure-drop distribution within the overall 
flow system. This reduction of orifice diameter simulates a change in manifold pressure 
losses. Such losses are frequently encountered in injector designs and contribute to sys- 
tem stability. 

Splash plate . - A splash plate was affixed to the reference injector (fig. 6(g)) to pro- 
duce a large change in the atomization and dispersion of the oxygen. The splash plate 
was used to test the effect on system stability of a significant change in the combustion 
process. 


EXPERIMENTAL EVALUATIONS 
Combustor Test Procedure 

Combustor performance was evaluated during runs in which propellant mass flow 
rates, supply pressures, and static as well as dynamic pressures were the primary per- 
formance measurements. The procedure for obtaining data was as follows: Hydrogen 
flow rate was held constant and oxygen flow rate was varied until the region between 
stable and unstable operation was well defined. This procedure was repeated at four dif- 
ferent hydrogen flow rates until a stability boundary was established for each configura- 
tion. 

Overall flow- system pressure drops were measured for each operating condition. 
Injector pressure drops used in this study were established in nonburning tests, and the 
deviation from burning conditions is unknown. Injector pressure-drop measurements 
attempted during test firings were unreliable because (1) the low pressure drops of the 
injector could not be established accurately at the operating pressure of the system and 
(2) the insertion of a transducer in the manifold caused changes in the stability of the 
system. 

Liquid-oxygen temperature was measured during burning tests by inserting a thermo- 
couple in the injector manifold. The average temperature was -305° F, which corre- 
sponds to a density of 0. 0421 pound per cubic inch. 


Experimental Results 

Experimental data for the reference combustor configuration are shown in figures 7 
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to 9. Combustion pressure is displayed in figure 7 as a function of the ratio of mass flow 
rates o/f for the four hydrogen flow rates. Pressure drops across the hydrogen injec- 
tion system exclusive of the critical-flow nozzle are also shown. A stability boundary 
was established by a line separating the stable and unstable conditions into two regions. 
Figure 8 is a cross plot of these results where combustion pressure and the stability 
boundary are presented as functions of oxygen and hydrogen flow rates. Also shown on 
this performance map are the oxygen injector-orifice pressure drop AP^ total oxygen 
flow-system pressure drop AP^., and the oscillatory frequency of the system at conditions 
along the stability boundary. Figure 9 shows the experimental data used to establish 
these oscillatory frequencies. Figure 10 is a performance map that displays computed 
values of c* efficiency, r Q c> over the range of test conditions. 

Experimental data for all combustor configurations were evaluated in a similar man- 
ner. The performance maps for these configurations are shown in figures 11 and 12. 


COMPARISON OF EXPERIMENTAL RESULTS 
Monopropellant Behavior 

In comparing the experimental results of the various configurations, it will be as- 
sumed that perturbations in hydrogen flow rate due to combustor pressure oscillations 
have a negligible effect on stability compared with the effect of perturbations in oxygen 
flow rate. Relative insensitivity to hydrogen flow perturbations can be justified for the 
injector configurations and mixture ratios used in these tests. Gaseous hydrogen was 
diffusely injected upstream of the oxygen orifices at a relatively high pressure drop and 
at a rate exceeding stoichiometric proportions. With such gaseous injection, any per- 
turbation in flow rate at the point of injection diffuses rapidly and does not at some later 
time cause a corresponding flow perturbation at the position of burning; however, a 
small perturbation in the total mass of hydrogen in the combustor can occur. Such a per- 
turbation has a reduced and negligible effect on combustor pressure because hydrogen in 
excess of that needed for chemical conversion prevails at all times. The dynamic behav- 
ior of the combustor system, therefore, is predominantly that of a monopropellant (oxygen) 
system with the bipropellant effects confined to establishing mean values of pressure, 
flow rate, and other parameters. 


Method of Comparison 

Although various comparisons of the performance maps of figures 11 and 12 can be 
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made, the experimental effect of configuration changes on stability will be examined by 
using properties that characterize the oxygen flow system and the combustion system. 
The properties and the reason for comparing these properties are given in the following 
paragraphs. 

Oxygen flow rates . - One measure of the stability of a combustor is the range of op- 
erating conditions that provides stable operation. The oxygen flow rate at the stability 
boundary is an index of the extent of the stable operating region; that is, a configuration 
change causing a decrease in this oxygen flow rate indicates an increased range of stable 
operation. For simplicity, the comparisons are restricted to the oxygen flow rate at a 
combustor pressure of 350 pounds per square inch absolute along the stability boundary; 
this is a condition that adequately represents each configuration. 

Mixture ratio . - Mixture ratio is a generally used performance parameter and is re- 
lated to oxygen flow rate at a pressure of 350 pounds per square inch absolute along the 
stability boundary. This parameter will also be used to indicate the extent of the stable 
operating region. As in the case of flow rate, a decrease in mixture ratio indicates an 
improvement in stability. 

Total pressure-drop parameter . - For steady flow, the oxygen flow process can be 
characterized by the relation 


AP, 



= K m 


which for small perturbations gives the relation 


m o - m O 


m. 


2 AP. 



The parameter P c /(2 AP t ) is an index of the magnitude of a flow perturbation caused 
by a pressure perturbation and, when sufficiently small, will generally ensure stability 
in a combustor system. The values of P c /(2 AP t ) required for stability will differ for 
various combustor systems because of interacting effects that are inherent in the overall 
system. A comparison of P / (2 A P t ) at the stability boundary of two combustors, there- 
fore, indicates the difference in their inherent stability; the combustor with a large value 
of P c /(2 AP t ) is inherently more stable. The configurations are compared by using this 
index, with the comparison again being restricted to a combustor pressure of 350 pounds 
per square inch absolute. 

Resistance parameter . - The resistance parameter is a measure of the total 
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resistance to steady flow of the oxygen flow system. Changes in this parameter indicate 
changes in the flow system. Although is a function of the flow area, its value is af- 

fected also by pressure recovery factors, coefficients of discharge, and viscous and other 
flow losses. Values of for different configurations will be compared to indicate 

changes in resistance; a large value of K m indicates a large resistance. 

Coefficient of discharge . - Values of the coefficient of discharge of the injector ori- 
fices C Di will also be compared. Although C^ is incorporated in K m , a change in 
C Di indicates a change in the discharge of oxygen from the injector that could affect pro- 
pellant preparation and the combustion rate. 

Combustor efficiencies . - Combustor efficiency is an index of the burning process. 

In this study, two variables that have important effects on burning rate are (1) a change 
in configuration and (2) a change in propellant flow rate. The efficiency if evaluated at 
mean oxygen and hydrogen flow rates of 0. 3 and 0. 15 pound per second, respectively, 
are compared to indicate the effect of a configuration change on the burning process inde- 
pendent of changes in stability. The efficiency 77 c evaluated at a combustor pressure of 
350 pounds per square inch absolute at the stability boundary will also be considered, be- 
cause flow rates at this condition vary substantially with changes in stability. This effi- 
ciency is a measure of the burning process for the flow-rate conditions that prevail at the 
stability boundary. 


Effect of Combustor Configuration Variables on Stability 


In the following tables, the previously discussed parameters for each configuration 
variable are compared with those for the reference combustor configuration. An in- 
crease or decrease in stability relative to the reference configuration is tabulated. The 
range of .the stable operating region and the inherent stability, indicated by m^ and 
P c /2 APj., respectively, are used in the stability rating. The tables are generally self- 
explanatory; however, several remarks about flow and combustion- system behavior are 
included. 

Characteristic length of combustor: 


Configuration 


System parameter 



Stability 

variable. 










L* 

m o 

o/f 

P c /(2 AP t ) 


C Di 

’'c 


Range 

Inherent 

81 

0.417 

3.35 

1.79 

563 

0. 505 

77 

78 

Decrease 

Decrease 

156 (reference) 

.37 

3.47 

2.24 

563 

. 505 

85 

88 



224 

.36 

3. 18 

2.34 

563 

. 505 

87 

90 

Increase 

Increase 
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The combustion process changed with L* because of the change in the length and volume 
for combustion. 

Length to diameter ratio of injector orifice: 


Configuration 


System parameter 



Stability 

variable, 

L/D 

“o 

o/f 

P c /(2 AP t ) 


C Di 



Range 

Inherent 

1. 56 (reference) 

0. 37 

3.47 

2.24 

563 

0. 505 

85 

88 





3.90 

.33 

2.65 

2. 57 

570 

.620 

92 

93 

Increase 

Increase 

6.25 

.335 

2.65 

1 

2.82 

552 

.675 

88 

89 

Increase 

Increase 


The combustion-process and injector-discharge changes with L/D may be related effects. 
Number of jets: 


Configuration 

System parameter 

Stability 

variable, 
number of jets 

m 0 

o/f 

P c /(2 AP t ) 

K 

m 

C Di 



Range 

Inhererit 

20 

40 (reference) 
60 

0.33 

.37 

.415 

2.36 

3.47 

4.80 

2.24 

2.24 

1.67 

715 

563 

610 

0. 509 
. 505 
. 568 

87 

85 

84 

87 

88 
i 80 

i 

Increase 

No change 

Decrease 

Decrease 


The flow-system resistance increased with 20 jets because of the reduced flow area. 
Injector-discharge, combustion- process, and flow-system-resistance changes with 60 
jets indicate a significant change in flow- system behavior and propellant preparation. 

Flow- system volume: 


Configuration variable, 


System parameter 



Stability 

flow system volume 

m 0 

o/f 

P c /(2 AP t ) 

K m 

C Di 

*>0 

*>0 

Range 

Inherent 

Reference 

1 0. 37 

3.47 

2.24 

563 

0. 505 

85 

88 





Large injector cavity 

.315 

2.3 

2.36 

745 

.680 

90 

91 

Increase 

Increase 

Extended manifold 

.38 

i 

3.6 

1.77 

685 

. 505 

83 

82 

Decrease 

Decrease 


The increase in flow-system resistance with an increase in volume may indicate a change 
in gas entrainment due to heat transfer and/or cavitation. The change in the combustion 
process may be related but shows no consistent trend with the volume change. 
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Diameter of upstream orifice:. 




Configuration 
variable, 
diameter of 
upstream orifice 

System parameter 

Stability 

“o 

o/f 

P c /( 2 AP t ) 


C Di 

r 'c 


Range 

Inherent 

0.052 

.070 (reference) 

0.24 

.37 

0. 94 
3.47 

1.23 

2.24 

2460 

563 

0. 505 
. 505 

94 

85 

CD 00 
00 00 

Increase 

Decrease 




The flow- system resistance increased because of the small orifice. 
Splash plate: 


Configuration 
variable, 
splash plate 

System parameter 

Stability 

m Q 

o/f 

P c /( 2 AP t ) 

K 

m 

C Di 



Range 

Inherent 

Without splash plate 

0.37 

3. 47 

2.24 

563 

0. 505 

85 

88 



(reference) 










With splash plate 

.22 

.96 

3. 58 

102 

. 505 

98 

94 

Increase 

Increase 


The change in propellant preparation was effective in changing the combustion pro- 
cess. Flow- system resistance was affected by the splash plate. 


Composite Comparison of Combustor Configuration Variables 

In the following table, the configuration variables are listed in the order of their rel- 


Configuration variable 


With splash plate 
Length to diameter ratio, 6.25 
Length to diameter ratio, 3. 9 
Large injector cavity 
Characteristic length, 224 in. 

20 Jets 
Reference 

Characteristic length, 81 in. 

Extended manifold 
60 Jets 

Upstream orifice diameter, 0. 052 in. 


Inherent 

stability, 

P c /(2 AP t ) 

Oxidant to fuel | 
mass flow 
rate ratio, 
o/f 

Rank 

according 
to stability 
rating 

3. 58 

0. 96 

2 

2.82 

2. 59 

5 

2. 57 

2.73 

6 

2.36 

2.30 

3 

2.34 

3.14 

7 

2.24 

2.36 

4 

2.24 

3.47 

9 

1.79 

3.35 

8 

1.77 

3.60 

10 

1. 67 

4.80 

11 

1. 23 

0.92 

1 
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ative inherent stability. The range of the stable operating region (as indicated by the o/f) 
is shown, and the rank according to this stability rating is noted. The 0. 052 -inch 
upstream-orifice configuration is the only configuration that is rated significantly differ- 
ent by the two rating methods. This comparison of the configurations shows the design 
features that provide the greatest stability. However, the validity of combining all these 
design features in a single configuration to obtain stability must be questioned because in- 
teracting effects can occur. Since these experimental effects have limited application to 
other combustors, additional qualitative and quantitative evaluations were made by em- 
ploying an analytical model to interpret these data. 


THEORETICAL ANALYSIS OF RESULTS 

The existing analytical models for low-frequency combustion instability (refs. 1 to 6) 
are interrelations of response functions for both the flow and combustion systems. In 
these models, a combustion time lag is generally postulated to characterize a property of 
the combustion process. In reference 1, a more realistic representation is attained by 
also postulating an interaction index. This index accounts for variations in combustion 
rates caused by physical and chemical kinetic changes due to chamber and injector geom- 
etry, types of propellants, etc. , between one combustor design and another. In the pres- 
ent analysis, the analytical derivations presented in reference 1 are used. The analytical 
description of the flow system is expanded, however, to include the effect of a discrete 
feedline pressure drop, as represented by the upstream orifice in the experimental con- 
figurations. 


Analytical Model 

A schematic diagram of the oxygen flow system considered for analytical study is 
shown in figure 13. The propellant flow originates from a constant-pressure supply tank 
and incurs orifice-type pressure drops at both the tank outlet and the entrance to the com- 
bustion chamber. Liquid capacitance and inertia are considered to be important only in 
the flow line between the two orifices. However, the capacitance is assumed to exist at a 
point source. 

The dimensionless equation that describes the flow- system dynamics is given as 
follows: 
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&U + £sd~ + /<?Y-— j cp + Jl + jd + (gjd + /) — + [//(l - Y)^ + /*y] — - 
\ ^ dz / ^ dz 2 


+ / 2 <?Y(1 - y)A. 

dz J 


K = 0 (1) 


where S, /> and are the elasticity, inertia, and pressure-drop parameters, respec- 
tively, of the flow system and are defined by the following relations: 


2 AP.p£Ae 


m0 


g 


2 A PjgA0g 




2 AP, 


The resistance parameter sd, which is defined as 


sd 


_ / C P i A i\ 2 
^Du^u/ 


( 2 ) 


is the modifying factor between equation (1) and the equation presented in reference 1. 
Equation (1) relates the fractional perturbation in propellant injection p^ to the fractional 
perturbations in chamber pressure cp. 

The dimensionless equation that describes the combustion system dynamics as de- 
rived in reference 6 is 


+ cp = p! + n(cp - cp') (3) 

dz 

where n is the pressure index of interaction and the primed parameters are functions of 
a retarded time, z - 6. By following the procedure of reference 1, the magnitude of the 
perturbations p 4 and cp in equation (1) are assumed to be the same as the primed term 
perturbations in equation (3), and the two equations are combined by elimination of p!. 
Thus, 
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(l-n + ^ = 


f*>(l + */) A + ft Y — 
<* dz 2 


1 + sd + (6 si + /) 


+•[/*(! - YW+ /^Y] A_ 




dz 


dz 


+ / 2 £Y(1 - Y) 

dz 3 


><P' - A <P' 
n 


( 4 ) 


This equation is of the form 


L {ty - ^ <5) 

and from reference 6 has a solution in Laplace notation of the form 

L 1 (s) = e“ s6 L 2 (s) (6) 


Solving equation (4) by comparison with equation (6) yields the following general equa- 
tion which gives the dynamic relation between the flow system and the combustion system: 

(1 - n + s + ne _s6 ) jl + si + ijSst + /)s + [ /g{l - Y )ji + /<?y]s 2 + / 2 S Y(1 - Y)s 3 } 

+ ?(1 + Sds + /<?Ys 2 )e“ s6 = 0 (7) 

The condition of neutral stability is defined when the real part X of the complex con- 
jugate root s = X + ia> is equal to zero. Solutions for neutral stability are obtained when 
the real and imaginary parts of equation (7) are equated to zero. Two specific solutions 
are of interest in this study; one is for the capacitance concentrated at the injector 
(Y = 1), and the other is for the capacitance at the upstream orifice (Y = 0). For Y = 1, 


2 2 r 

(1 - n) 2 + co 2 Cs* + 1 - /<?w 2 ) +{sigu + /w) 2 J = |n^ + (n + ^)(l- /Ao 2 ) + n/co + (n+ 


6 = — <7r - tan 
Cl) 


• Y-"A - tan V y*- + /a. \ + tan - 1 
1 - n V+l-r#u> 2 / 


(n + P)si£oJ + n/to 


nst + (n +i?)(l - 


(8) 


( 9 ) 
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and for Y = 0, 



- fs4£i>)) + {sd£ co + 


/-) 2 ] 


8 =ll - tan" 1 \ - tanV dS(a + >— \ 

\1 - n/ \s4 + 1 - /j*gw 2 ) 

For each of the two cases, solutions for 6 
both corresponding equations. 


= n(l + J - Jrfgu 2 ) + <?j 2 

+ [(n + + n/w] 2 

n Ju \ 


+ tan 


(n + &)si£ <x> + 


( 10 ) 

(ID 


[n {sd + 1 - JsiS to ) + ^JJ 
are obtained by simultaneously satisfying 


Analytical Solutions 

Solutions for 5 along the stability boundaries shown in figure 11 are only possible 
for known values of either n or g . The interaction index n characterizes the combus- 
tion dynamics and is not specified by the experimental measurements. The elasticity g 
is considered unknown because of the probability of gas entrainment in the cryogenic flow 
system. Discrete points along the stability boundaries, therefore, have a range of prob- 
able solutions for 5 . Solutions for 6 and g for assumed values of n were obtained by 
iteration on a high-speed digital computer. The related dimensional quantities of time 
lag r and liquid compressibility e for these solutions are given in table I. 

For a typical operating condition figure 14 shows the relation between n, r, and e 
given by these solutions. Solutions for Y = 0 and Y = 1 are shown. For both these 
solutions, an upper- limit value of n is obtained when liquid compressibility is equal to 
zero. Realistically, n must be less than the value corresponding to liquid compressibil- 
ity for oxygen with no gas entrainment (10x10” psi ), as noted in figure 14. Other op- 
erating conditions shown in table I exhibit a similar relation between n and t; however, 
the limiting value of n for some conditions is less than zero. Negative values of n 
give valid solutions for these conditions, but the numerical evaluations were not made in 
this study. 


DISCUSSION OF ANALYTICAL SOLUTIONS 
Position of Flow-System Elasticity 

Table I presents solutions for the reference combustor configuration with the flow- 
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system elasticity assumed concentrated at the injector Y = 1 and at the upstream orifice 
Y = 0. Figure 14 shows these solutions for one operating condition. Differences in the 
solutions for Y = 1 and Y = 0 are negligible for this and all other configurations. It is 
concluded that stability is insensitive to the position of any local cavitation or gas entrain- 
ment in the flow system. The effective position of such a gas inclusion, therefore, cannot 
be identified from the stability data. 

This insensitivity to position is attributed to the dynamics of the particular flow sys- 
tem used for the experimental tests. An evaluation of the phase- angle terms in equa- 
tions (9) and (11) shows that the delay between a combustor pressure change an a flow 
change is small compared with the period of oscillation. This causes the stability of the 
overall combustor system to be relatively insensitive not only to the position but also to 
small changes in flow- system capacitance. 


Composite Evaluation of Interaction Index and Time Lag 

Since both the flow- system elasticity and the interaction index are unknown, unique 
values of t, e, and n for each test condition cannot be established. Some insight into 
the behavior of the combustion and flow systems can be acquired, however, by analyzing 
the solutions for two conditions; these are (1) a constant liquid compressibility and (2) a 
constant interacting index. 

Interaction index . - Values of interaction index are obtained for each operating con- 
dition when liquid compressibility is taken as that of liquid oxygen without gas entrain- 
ment (10x10” psi ). These values of interaction index range between approximately 
0. 7 and 0 and are shown as a function of combustor pressure in figure 15. Solutions for 
the splash-plate configuration were not evaluated but are known to be negative. The 
curves in figure 15 represent individual configurations. Although some increase in the 
index with pressure is noted for each configuration, the interaction index does not appear 
to be a unique function of pressure. The dispersion of values shown in figure 15 implies 
that the pressure sensitivity of the combustion process as indicated by the parameter n 
must vary with operating conditions and configuration changes in order that the assump- 
tion of a constant compressibility be valid. 

A similar dispersion of results is obtained for a constant compressibility other than 
that for liquid oxygen. Assuming a liquid compressibility greater than that of liquid oxy- 
gen gives smaller values of interaction index, as shown in table I and demonstrated in fig- 
ure 14. The values of interaction index shown in figure 15, therefore, are the highest 
possible values for each condition (since compressibility less than that for liquid oxygen 
is impossible). Hence, any comparison of results made on the basis of a constant inter- 
action index must be made for an index at or near zero if the comparison is to encompass 
all test conditions. 
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Probable variations in liquid compressibility may be examined in a similar manner. 
Values of liquid compressibility are obtained for each operating condition when interac- 
tion index is assumed constant. Figure 16 shows these values of liquid compressibility 
for n = 0 as a function of injector manifold pressure with the curves again relating data 
for each configuration. Although a large dispersion of results is obtained, a decrease in 
compressibility with an increase in pressure is apparent when characteristics of individ- 
ual configurations are analyzed. This trend is consistent with the properties of a two- 
phase medium, as shown by the lines of constant volume percent gas in the flow system 
in figure 16. Gas entrainment of less than 2 percent of the flow- system volume could ac- 
count for the values of compressibility needed to satisfy analytical solutions for an inter- 
action index of zero. 

Gas entrainment in the flow system was confirmed in several tests. The injector 
cavity was observed through window ports, and a dispersion of small bubbles was noted. 
Similar behavior in a chugging combustor was noted in reference 7. Although quantitative 
measurements from these observations were not possible, it may be concluded that gas 
entrainment modified the elasticity of the flow system. Therefore, the expected values 
of interaction index are less than those shown in figure 16, which assume no gas entrain- 
ment. The obvious conclusion is that the interaction index for these combustors exhibits 
a dispersion of values about a mean level of less than 0. 25 and probably near zero. 

Time lag . - Time lag is relatively insensitive to the specific parametric values used 
to satisfy analytical solutions for each test condition. This insensitivity is shown in fig- 
ure 17, where time lag is presented as a function of combustor pressure for the condition 
of interaction index equal to zero and also for compressibility equal to that for liquid oxy- 
gen. Values are somewhat larger when the index is zero. Although a decrease in time 
lag with an increase in pressure appears to be a characteristic property of each configur- 
ation, the dispersion about a mean variation with pressure is relatively large. Appar- 
ently, a unique time lag and pressure relation cannot be used to characterize these con- 
figurations. 

A variation in time lag with several parameters may be expected. The time lag rep- 
resents the average time interval between propellant injection and complete chemical con- 
version. The time required for propellant atomization is included in this time interval. 

In this study with parallel-jet injection methods, the jet-breakup, or atomization, time 
is a substantial portion of the total time lag. Parameters affecting jet-breakup time, 
therefore, must be included in any correlation of time-lag variations. 

Although a precise value of time lag cannot be established for each test condition, 
the values shown in figure 17 exhibit trends that may be related to changes in jet- breakup 
time. A degree of consistency exists when time lag is assumed to be equal to the summa- 
tion of a jet-breakup time, vaporization time, and a gas- phase- mixing time. When the 
mean vaporization and mixing times are approximated by the theoretical combustion gas 
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residence time 0 g , jet-breakup time is given by r - 0 g . In the subsequent discussion, 
changes in the value of r - 0 g are compared with the characteristics of a jet-breakup 
process. The comparison shows a relation between changes in stability and changes in 
one mechanism within the overall combustion process. 


EVALUATION OF JET-BREAKUP TIME 
Breakup Processes 

Jet breakup has been studied by various investigators (refs. 8 to 10) who showed that 
breakup in the concurrent flow of a gas and a liquid is related to both the velocity of the 
liquid jet and the relative velocity between gas and liquid. In general, breakup time in 
the turbulent flow region decreases with an increase in jet velocity and in relative veloc- 
ity, where the increased relative velocity is larger at high gas pressures. The effect of 
jet velocity on breakup time varies with changes in the jet-formation process. The jet 
may fill the exit orifice, partially separate from the sidewalls of the orifice, or flow com- 
pletely free of the sidewalls. Partial separation is analogous to cavitation collapse and 
reduces breakup time relative to breakup when the jet is flowing freely or full. Partial 
separation and flowing full exhibit higher flow coefficients than that exhibited by flowing 
free because of changes in effective flow area. Increasing the orifice length of a free- 
flowing configuration initially causes a transition to partial separation and finally to flow- 
ing full. Changes in the liquid cross velocity at the orifice inlet affect the conditions at 
which these transitions occur. 

These qualitative properties of jet breakup will be used to interpret changes in time 
lag caused by probable changes in jet- breakup time. By this interpretation, the stability 
characteristics of the various configurations will be related to the jet-breakup process, 
where breakup is assumed to be controlled by flow conditions prior to any significant 
burning of propellants. 


Effect of Combustor Configuration Variables on Breakup Time 

The properties of each configuration at a combustor pressure of 350 pounds per 
square inch absolute along the stability boundary are tabulated and compared with those 
for the reference combustor. Parameters used to characterize jet breakup are the oxy- 
gen jet velocity Uq, the injector orifice discharge coefficient C Di , and the relative ve- 
locity term p g /g(U H - Uq) 2 expressed as a dynamic pressure. Oxygen jet velocity is 
the computed value based on mass flow rate and orifice exit area. Hydrogen velocity is 
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mu 


in ii 


i ii 


the concurrent flow velocity surrounding the oxygen jet, computed from the mass flow 
rate and from the combustor cross-sectional area at inlet temperature conditions. The 
term t - is used to indicate the apparent jet-breakup time, as previously discussed. 
Time-lag values obtained from solutions for an interaction index of zero are used, al- 
though the comparison is not dependent on the values of this specific solution. 

Characteristic length of combustor: 


Configuration 

variable, 

L* 

System parameter 

Jet- 

breakup 

time 


u o 

C Di 

U H 

(p g /g)(U H - U Q ) 2 

0 g 

T 

T -* g 

81 

16. 5 

0. 505 

52.8 

4.45 

0.0025 

0. 0189 

0.0164 

156 (reference) 

14.7 

. 505 

45. 5 

3.21 

.0044 

.0228 

.0184 

224 

14.2 

. 505 

47.3 

3. 50 

.0061 

.0243 

.0182 


The changes in jet-breakup time r - 6 ^, although small, show a decrease with an in- 
crease in jet velocity and dynamic pressure; this behavior is expected for a jet-breakup 
process. However, the change in vaporization and mixing times 9 appears to cause a 
larger change in time lag and presumably in stability than that caused by the change in jet- 
breakup time. The constant discharge coefficient Cp. implies that no change in jet 
formation process is encountered. 

Length to diameter ratio of injector orifice: 


Configuration 

variable, 

L/D 



System parameter 



Jet- 

breakup 

time 

u o 

C Di 

U H 

V g)(U H' U 0 )2 

*g 

T 

T -*g 

1. 56 (reference) 

14.7 

0. 505 

45. 5 

3.21 

0. 0044 

0. 0228 

0.0184 

3.90 

13. 1 

.620 

53. 1 

5.40 

.0039 

.0185 

.0146 

6.25 

13.2 

.675 

53.6 

5. 51 

.0042 

.0185 

.0143 


An increase in orifice L/D decreased the jet-breakup time. The decrease is not 
consistent with the behavior expected for the decrease in jet velocity, or for the increase 
in dynamic pressure. A change in jet formation, however, is probable, as indicated by 
the increase in the discharge coefficient. This increase is characteristic of a jet that 
flows freely at a small L/D but is only partially separated from the orifice sidewalls at 
a larger L/D. The more turbulent jet formation with partial separation is a probable 
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cause of the decrease in breakup time. An examination of all the data showed that this 
turbulent jet was less sensitive to changes in dynamic pressure than a free jet. 

Number of jets: 


Configuration 
variable, 
number of jets 



System parameter 



Jet- 

breakup 

time 


u o 

C Di 

U H 

(Pg/g )(U H-U 0 ) 2 

0 g 

T 

T - 0 g 

20 

1 

26.2 

0. 509 

59 

3.64 

0. 0042 

0.0181 

0.0139 

40 (reference) 

14.7 

. 505 

45. 5 

3.21 

.0044 

.0228 

.0184 

60 

11. 1 

. 568 

37. 1 

2.30 

.0048 

.0191 

.0143 


Jet-breakup time is smaller for both the 20- and the 60-jet configurations than for 
the 40-jet reference configuration. The decrease in breakup time for the 20-jet config- 
uration appears consistent with the increase in jet velocity. The decrease for the 60-jet 
configuration, however, must be attributed to a change in jet formation, as implied by the 
increase in the discharge coefficient. The relative cross velocity is higher for the 60-jet 
configuration than for a configuration with fewer jets and may have promoted partial flow 
separation from the orifice sidewalls. 

Flow- system volume: 


Configuration variable, 
flow- system volume 

System parameter 

Jet- 

breakup 

time 

u o 

C Di 

U H 

(Pg/g )(u H -U 0 ) 2 

0 g 

T 

T -*g 

Reference 

14.7 

0. 505 

45. 5 

3.21 

0. 0044 

0.0228 

0.0184 

Large injector cavity 

12. 5 

.680 

58.2 

7.06 

.0040 

.0181 

.0141 

Extended manifold 

15.1 

.505 

44.7 

2.96 

.0050 

.0242 

.0192 


The small increase in breakup time for the extended manifold configuration is in 
agreement with the small changes in the velocity and dynamic pressure parameters. The 
decrease in breakup time for the large injector cavity configuration is significant; it may 
be partially attributed to the large increase in dynamic pressure. However, the increase 
in discharge coefficient indicates a condition of partial flow separation, which would also 
decrease breakup time. Lower cross velocities with the large cavity could cause such a 
jet formation change. 
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Diameter of upstream orifice: 


Configuration variable, 
diameter of upstream 
orifice 



System parameter 



Jet- 

breakup 

time 


u o 

C Di 

U H 

(p g /g )(u H - u 0 ) 2 

e g 

T 

T - e g 

0.052 

9. 5 

0. 505 

110.0 

34.2 

0. 0040 

0.0150 

0.0110 

. 070 (reference) 

14.7 

.505 

45. 5 

3.21 

.0044 

.0228 

.0184 


Flow conditions at the stability boundary were significantly different for the small up- 
stream orifice from those for the reference configuration. The order- of-magnitude in- 
crease in dynamic pressure apparently reduced the jet-breakup time even though the de- 
crease in jet velocity would cause the opposite effect. Jet-formation process appears un- 
affected. 

Splash plate: 

A numerical comparison between the splash- plate configuration and the reference 
configurations is not possible because the solutions for a negative interaction index re- 
quired for the splash plate were not evaluated. 


Effect of Jet Breakup on Interaction Index 

The interaction index n is a measure of the pressure sensitivity of the combustion 
process; it is also a measure of the perturbation in time lag caused by a perturbation in 
combustor pressure. It has been previously established that a value of n at or near 
zero is the probable average value for these configurations. Pressure oscillations, 
therefore, should cause small variations in jet-breakup time, since delay time and break- 
up time are closely related. The qualitative properties of jet breakup do not clearly indi- 
cate a small value of n. However, an increase in combustor pressure would decrease 
jet velocity and increase delay time (negative values of n), whereas an increase in pres- 
sure would also increase the dynamic pressure and decrease breakup time (positive 
values of n). The relative importance of these two effects would require a quantitative 
description of the jet-breakup process. 

In the case of the splash-plate configuration, where a negative interaction index is 
clearly established, a decrease in time lag with pressure should be evident. Here, jet 
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impingement is on a solid body, and breakup or atomization time may be more dependent 
on jet velocity than it is for the nonimpinging jets. This would reduce the effectiveness 
of the dynamic pressure on atomization and cause the negative n characteristic of a jet- 
velocity increase to predominate. 


AVERAGE FUNCTIONAL VARIATION OF TIME LAG 


The jet formation of the reference configuration is believed to be that of a jet flowing 
freely of the orifice sidewalls. This jet form was exhibited by several configurations that 
provided stability data over a broad range of operating conditions. These data were used 
to formulate a relation for jet breakup. 

In formulating the relation, a basic breakup time was assumed to be established by 
jet turbulence and other flow disturbances caused by flow through an orifice. This basic 
breakup time was reduced by the dynamic pressure forces acting on the jet during the 
period of breakup. The relation formulated on the basis of such a model is 


T i = 


0.2 


U 


0.75 


- T i 


Per 9 

1 T (U H - u o> 


or 



U 


0.2 

0.75 

O 


1 + 0 . 


p 

1-S (U : 


H 


u o) 


2 


With the aid of this equation for jet-breakup time and the relation r = Tj + 0 g , 
values of time lag are defined for all operating conditions of the reference combustor 
configuration. Figure 18 shows time lags for the range of conditions considered in this 
study. In general, time lag varies with hydrogen flow rate and is relatively constant with 
variations in oxygen flow rate. 

Figure 18 displays the time lags that characterize the combustion process for a 
range of operational conditions. The specific values that apply to the condition of neutral 
stability depend on the interaction of many stability parameters and analytically involve a 
complex interactive procedure where boundary conditions continually change until conver- 
gence at neutral stability is obtained. Such solutions are possible; however, greater 
precision is needed in describing time-lag variations with operating conditions if these 


22 



variations are to be useful in analytically modeling the effect of system variables on sta- 
bility observed in this study. 


SUMMARY OF RESULTS 

Chugging instability limits were experimentally established in a parametric study of 
a small-scale combustor system that used parallel- jet injection of liquid oxygen into dif- 
fusely injected gaseous hydrogen. System stability was experimentally characterized by 
the range of the stable operating region and by the inherent stability related to the oxygen 
flow- system drop at the stability boundary. An analytical model of the experimental com- 
bustor system was developed. With this model the configurations were characterized by 
a time lag (delay time between injection and burning) and by an interaction index (sensitiv- 
ity of combustion to pressure oscillations). 

In general, no predictable trends in the variations of interaction index were noted. 
The average value for all tests was less than 0. 25. Time lag appeared to vary systemat- 
ically and was related to the sum of a liquid jet-breakup time and a vaporization and mix- 
ing time, with the contribution of breakup time being predominant. Breakup time was de- 
pendent on the jet-formation process (flowing freely or attached to the orifice sidewalls) 
and on the aerodynamic forces imposed by the combustor gases. A relation for breakup 
time was formulated for a jet flowing freely of the sidewalls. 

The principal effects of specific combustor-system variables on stability are sum- 
marized as follows: 

1. Characteristic length: An increase in characteristic length increased the range of 
stable operation, inherent stability, time lag (vaporization and mixing), and combustor 
performance. Jet flow remained free of the orifice sidewalls. 

2. Injector orifice length to diameter ratio: An increase in orifice length to diameter 
ratio increased the range of stable operation and the inherent stability. Jet flow appar- 
ently attached to the sidewalls with an increase in length to diameter ratio causing a de- 
crease in time lag (breakup time) and an increase in flow- system pressure drop. 

3. Number of jets: An increase in the number of injector orifices, or jets, de- 
creased the range of stable operation and inherent stability. A large number of jets af- 
fected the jet-formation process and decreased the time lag (jet-breakup time). 

4. Flow-system volume: Increasing the flow-system volume by extending the injec- 
tor manifold decreased the range of stable operation and the inherent stability, whereas 
increasing the size of the injector cavity had the inverse effect. Gas entrainment due to 
heat transfer and/or local cavitation was experimentally and analytically identified within 
the combustor system and may have contributed to these effects. Time lag (jet-breakup 
time) decreased with enlargement of the injector cavity, and jet flow exhibited a condition 
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of partial separation from the orifice sidewalls. 

5. Upstream orifice: Increasing the flow-system pressure drop by inserting a small 
orifice in the supply line increased the range of stable operation but decreased the inher- 
ent stability. Time lag (jet-breakup time) decreased in a predictable manner for a free- 
flowing jet. 

6. Splash plate: Impingement of oxygen jets on a splash plate increased the range 
of stable operation and the inherent stability. Negative values of the interaction index 
were established for the splash-plate configuration; however, time-lag values were not 
evaluated analytically. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 17, 1966, 
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TABLE I. - ANALYTICAL SOLUTIONS 


Operating parameter 

Interaction index, n 

Mass flow rate, 
lb/sec 

Com- 

bus- 

tor 

pres- 

sure, 

P c’ 

psi 

Charac- 

teristic 

exhaust 

velocity, 

c*, 

ft/sec 

Pressure - 
drop param- 
eter of in- 
jector 
orifice, 

9 

Fre- 

quen- 

cy, 

f, 

cps 

Mean 

gas 

resi- 

dence 

time 

in 

com- 

bus- 

tor, 

V 

msec 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

Liquid 

compress- 

ibility, 

€, 

psi" 1 

Time 

lag, 

T, 

msec 

Liquid 

compress- 

ibility, 

€, 

psi" 1 

Time 

lag, 

T, 

msec 

Liquid 

compress- 

ibility, 

e, 

psi -1 

Time 

lag, 

T, 

msec 

Liquid 

compress- 

ibility, 

> 

Time 

lag, 

T, 

msec 

Liquid 

compress- 

ibility, 

€, 

psi" 1 

Time 

lag, 

T, 

msec 

Liquid 

compress- 

ibility, 

€, 

psi" 1 

Time 

lag, 

T, 

msec 

Liquid 

compress- 

ibility, 

e, 

psi -1 

Time 

lag, 

msec 

Oxy- 

gen, 

m c 

Hydro- 

gen, 

m H 

Reference configuration 3 (Y = 0) 

0.323 

0.088 

288 

6770 

29.4 

20 

4.60 

62.9X10' 6 

28.1 

51.1X10" 6 

26.4 

35.8X10" 6 

24.2 

0. 815x10"® 

19.4 










.411 

. 129 

400 

7170 

25.6 

28 

4.30 

46.6 

19.9 

40.5 

18.5 

33.5 

17. 7 

24.4' 

16.2 

9. 50x10“® 

13.8 


— 


— 

.463 

.171 

487 

7440 

24.3 

34 

4.09 

38.0 

16.2 

33.9 

15.5 

29.2 

14.6 

23.6 

13.6 

16.4 

12.3 

0. 130X10" 6 

9.86 


— 

.497 

.207 

538 

7390 

23.6 

37 

4.06 

35.1 

14.8 

31.8 

14.2 

27.9 

13.4 

23.5 

12.6 

18.1 

11.5 

10.3 

10.1 


— 


Reference configuration 3, (Y = 1) 


0.323 

1 

0.088 

j 

288 

6770 

l 

29.4 

20 

4.60 

67.0X10" 6 

28.5 

54. 6X10"® 

26.8 

38.6X10" 6 

24.6 

5. 65x10"® 

20.0 









.411 

.129 

400 

7170 

25.6 

28 : 

4.30 

51.0 

20.3 

44.4 

19.3 

36.7 

18.1 

27.3 

16.5 

12.1X10" 6 

14.2 


— 


— 

.463 

.171 

487 

7440 

24.3 

34 

4.09 

42.2 

16.6 

37.7 

15.9 

32.6 

15.0 

26.5 

13.9 

18.8 

12.6 

4. 16X10" 6 

10.2 


— 

.497 

.207 

538 

7390 

23.6 

37 

4.06 

39.6 

15.2 

35.8 

14.5 

31.3 

13.8 

26.5 

12.9 

20.7 

11.8 

12.5 

10.4 


— 


Characteristic length, L*, 81 in. (Y = 0) 


0.354 i 
.422 

0.088 

.129 

255 

359 

5580 ! 

6310 

21.6 

21.6 

24 

34 

2.92 

2.54 

64. 5X10" 6 
39.6 

24.9 

17.2 

54.2x10"® 

33.8 

23.6 

16.3 

41.2X10" 6 

26.9 

22.0 

15.3 

21.8X10" 6 
17. 3 

19.5 

13.8 




— 


— - 

.469 

.172 

436 

6580 

21.4 

40 

2.40 

31.5 

14.4 

27.2 

13.8 

22.2 

12.9 

15.6 

11.7 

2.53X10" 6 

9.62 


— 


— 

.502 

.213 

494 

6680 

21.1 

46 

2.33 

26.6 

12.4 

23.3 

11.8 

19.4 

11.1 

14.6 

10.2 

7.30 

8.91 


... 




Characteristic length, L*, 224 in. (Y = 0) 


0.284 

0.088 

! 

: 275 

7180 
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Extended manifold (L = 70. 40, Y = 0) 
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a Length to diameter ratio of injector orifice, 1. 56; 40 jets; length of flow-system supply line, 22. 4 in. ; area of upstream orifice, 0. 00384 in. area of injector orifice, 0.0503 in. 2 ; 
characteristic length of combustor, 156 in. ; resistance parameter of flow system, 47. 348; effective area of flow system supply line, 0. 14 in. 2 . 

















(a) Reference configuration; 40 jets; length to diameter 
ratio, 1.56; cavity volume, 0. 25 cubic inch. 



(b) Length to diameter ratio, 3. 90. 



(c) Length to diameter ratio, 6. 25. 



(f) Cavity volume, 1. 56 cubic inches. 



(g) Splash plate. 

Figure 6. - Liquid-oxygen injector cavity and orifice configurations. 
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Figure 7. - Chamber pressure as function of oxidant-fuel ratio for reference combustor. 
40 Parallel jets; length to diameter ratio of injector orifice, 1.56; characteristic length 
of combustor, 156 inches; volume of liquid-oxygen cavity, 0. 25 cubic inch; diameter of 
upstream orifice, 0.070 inch. 
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Figure 8. - Composite stability plot of liquid-oxygen weight flow rate as function of gaseous-hydrogen weight 
flow rate for reference combustor. 40 Parallel jets; length to diameter ratio of injector orifice, 1.56; volume 
of liquid-oxygen cavity, 0. 25 cubic inch; characteristic length of combustor, 156 inches; diameter of up- 
stream orifice, 0.070 inch. 
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Figure 9. - Chugging frequency as function of 
chamber pressure. 40 Parallel jets-, length 
to diameter ratio of injector orifice, L 56; 
characteristic length of combustor, 

156 inches; diameter of upstream orifice, 
0.070 inch; volume of liquid-oxygen cavity, 

0. 25 cubic inch. 



Mass flow rate of hydrogen, m H , lb/ sec 

Figure 10. - Reference combustor efficiency map. 40 Parallel jets; 
length to diameter ratio of injector orifice, 1.56; characteristic 
length of combustor, 156 inches. 
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(c) Characteristic length of combustor, 
??A inrhes 




Combustor pressure, psia 
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\ Denotes value of oscil- 
latory frequency in cps 



Mass flow rate of hydrogen, m H , Ib/sec 

(j) Diameter of upstream orifice, (k) Splash plate. 

0. 052 inch. 



performance maps. 
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(a) Reference combustor. 



(d) Length to diameter ratio of injector 
orifice, 3.90. 



(h) Large injector cavity. 
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(b) Characteristic length of combustor, 
81 inches. 
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ft) Characteristic length of combustor, 
224 inches. 


Characteristic velocity, percent 
of theoretical 

Stability boundary; stable above, 
unstable below 



(e) Length to diameter ratio of injector 
orifice, 6.25. 
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(k) Splash plate. 


Figure 12. - Combustor efficiency maps. 
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Figure 13. - Schematic diagram of oxygen flow system. 




Figure 14. - Range of probable analytical solutions for 
one operating condition. 
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Figure 15. - Relation between interaction index and 
combustor pressure for liquid compressibility of 
lOxlO -6 psi" . 




Figure 16. - Liquid compressibility for solutions with 
interaction index of zero. 
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Combustor pressure, P c , psia 


(b) Interaction index, zero. 


Figure 17. -Time lag solutions for constant 
compressibility and constant interaction 
index. 



Figure 18. - Time-lag variations for burning rate process con- 
trolled by typical jet breakup. 
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"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the at?nosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof” 

—National Aeronautics and Space Act of 1958 
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